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DESCRIPTION 

CONTROL SYSTEM OF INTERNAL COMBUSTION ENGINE 
TECHNICAL FIELD 

The present invention relates to a control system of 
an internal combustion engine. 
BACKGROUND ART 

Regarding an internal combustion engine provided 
with a plurality of cylinders, introducing air into an 
intake pipe extending from a throttle valve to an intake 
valve through the throttle valve in exactly a throttle 
valve air passage amount, and discharging air from the 
intake pipe through an intake valve in exactly a cylinder 
air filling amount to fill a cylinder at the time of an 
intake stroke, known in the art is a control system of an 
internal combustion engine wherein a formula obtained 
from the Law of the Conservation of Mass for the intake 
pipe and a state equation for the air in the intake pipe 
is used to calculate the cylinder air filling amount of 
the cylinder (see Japanese Patent Publication (A) No. 
2002-70633) . 

To use the above formula to calculate a cylinder air 
filling amount, for example, the temperature of the air 
in the intake pipe and the volume of the intake pipe have 
to be found. However, to find the air temperature, for 
example, not only is a temperature sensor necessary, but 
also, if considering the response delay, even if using a 
temperature sensor, it would be difficult to accurately 
find the air temperature. Further, the intake pipe 
includes manufacturing error, so the volume of the intake 
pipe cannot be considered equivalent to, for example, the 
design value. Measurement of the volume for each intake 
pipe is also extremely impractical. 

Further, when using a formula to estimate a cylinder 
air filling amount, if using a formula obtained based on 
the Law of the Conservation of Mass and state equation as 



it is, the formula would become complicated and the 
calculation load would become huge, so normally this 
formula is simplified for use. Here, if the intake valve 
opening/closing timing is set to the retarded side, an 
intake valve will remain open even after intake bottom 
dead center. In this case, even if the piston starts to 
rise, since the intake valve is held in the open state, 
the air sucked into the cylinder is liable to flow back 
into the intake pipe. However, if simplifying the formula 
for use when estimating the cylinder air filling amount 
as explained above, this backflow of air is not 
considered and therefore the calculated cylinder air 
filling amount ends up including error. 
DISCLOSURE OF THE INVENTION 

Therefore, the present invention has as its object 
to provide a control system of an internal combustion 
engine able to simply and accurately calculate a cylinder 
air filling amount. 

The present invention provides as means for solving 
this problem a control system of an internal combustion 
engine as set forth in the claims. 

To solve this problem, in a first aspect of the 
invention, there is provided a control system for an 
internal combustion engine provided with a plurality of 
cylinders, introducing air into an intake passage part 
extending from a throttle valve to an intake valve 
through the throttle valve in exactly a throttle valve 
air passage amount, and discharging air from the intake 
passage part through an intake valve in exactly a 
cylinder air filling amount to fill a cylinder at the 
time of an intake stroke, wherein the cylinder air 
filling amount is divided into a first amount of air and 
a second amount of air, the first amount of air being an 
excess of a cylinder air filling amount with respect to a 
throttle valve air passage amount occurring due to an 
intake stroke, wherein the control system comprises an 
intake pressure drop detecting means for detecting a drop 



- 3 - 



in intake pressure occurring due to an intake stroke 
being performed for each cylinder; a first air amount 
calculating means for calculating the first amount of air 
for a cylinder based on its intake pressure drop; a 
5 throttle valve air passage amount detecting means for 

detecting a throttle valve air passage amount; a second 
air amount calculating means for calculating the second 
amount of air for a cylinder based on the throttle valve 
air passage amount; a cylinder air filling amount 

10 calculating means for totaling the first amount of air 

and the second amount of air to calculate the cylinder 
air filling amount for a cylinder; and a control means 
for controlling the engine based on the cylinder air 
filling amount of the cylinder, and wherein the first air 

15 amount calculating means sets a set crank angle range so 

as to include the intake strokes of at least two 
cylinders for which cylinder air filling amounts are to 
be calculated, calculates the total value of the intake 
pressure drop of the cylinders performing an intake 

20 stroke in the set crank angle range, and calculates the 

first amount of air based on each intake pressure drop 
and the intake pressure drop total value. 

Further, in a second aspect of the invention, there 
is provided the first aspect of the invention wherein 

25 when backflow of air from inside a cylinder to the intake 

passage part occurs at the end of an intake stroke, the 
action of the second air amount calculating means 
calculating the second amount of air is prohibited. 

To solve the above problem, in a third aspect of the 

30 invention, there is provided a control system for an 

internal combustion engine provided with a plurality of 
cylinders and a plurality of intake valves, wherein the 
cylinder air filling amount to a cylinder is divided into 
a basic amount of air and an excess amount of air flowing 

35 from an intake passage part to the inside of the cylinder 

exceeding a throttle valve air passage flow rate due to 
opening of an intake valve, and wherein the control 



system comprises a basic air amount calculating means for 
calculating a basic air amount based on a throttle valve 
air passage flow rate of air flowing into the intake 
passage part through the throttle valve and the opening 
time of each intake valve; an excess air amount 
calculating means for calculating an excess air amount 
based on the drop in intake pressure due to opening of 
the intake valve; a cylinder air filling amount 
calculating means for totaling the basic air amount and 
excess air amount to calculate a cylinder air filling 
amount to a cylinder; and a control means for controlling 
the engine, based on the cylinder air filling amount to a 
cylinder, and wherein the basic air amount calculating 
means calculates a virtual intake valve opening time so 
that the average air flow rate to all cylinders becomes 
equal to the throttle air passage flow rate and uses the 
virtual intake valve opening time as the opening time of 
an intake valve. 

According to the third aspect of the invention, the 
virtual intake valve opening time becomes a value by 
which the average air flow rate to all cylinders becomes 
equal to the throttle air passage flow rate. For this 
reason, when there is backflow of air from inside a 
cylinder to the intake pipe, the virtual intake valve 
opening time becomes shorter than the actual intake valve 
opening time. If this virtual intake valve opening time 
is used by the basic air amount calculating means to 
calculate the basic air amount, the basic air amount can 
be accurately calculated. 

Further, in a fourth aspect of the invention, there 
is provided the third aspect of the invention wherein the 
basic air amount calculating means uses the virtual 
intake valve opening time as the opening times of the 
intake valves when backflow of air to the intake passage 
part occurs near the intake valve opening timing or near 
the intake valve closing timing. 

According to the present invention, the cylinder air 
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filling amount can be simply and accurately calculated. 

The present invention will be more fully understood 
from the drawings and the description of preferred 
embodiments of the present invention. 
5 BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is an overall view of an internal combustion 
engine , 

FIG. 2 is a view of an intake valve opening/closing 
timing. 

10 FIG. 3 is a view of the results of detection of an 

intake pressure Pm. 

FIG. 4 is a time chart for explaining an intake 

pressure drop APmdwni. 

FIG. 5 is a view for explaining a method of 
15 calculation of a cylinder air filling amount Mci . 

FIG. 6A and FIG. 6B are time charts for explaining a 
method for calculation of a parameter Km. 

FIG. 7A and FIG. 73 are time charts for explaining 
another example of a method for calculation of a 
20 parameter Km. 

FIG- 8 is a time chart for explaining a cylinder 
intake air flow rate mci when an intake valve 
opening/closing timing is set to a retarded side. 

FIG. 9 is a flow chart showing a routine for 
25 calculation of a fuel injection time TAUi. 

FIG. 10 is a flow chart showing a routine for 
calculation of a cylinder air filling amount Mci for a 
first embodiment. 

FIG. 11 is a time chart for explaining error in 
30 approximation. 

FIG. 12A and FIG. 123 are time charts for explaining 
a method of calculation of a virtual intake valve opening 
t ime X . 

FIG. 13 is a time chart for explaining a method of 
35 calculation of the virtual intake valve opening time x. 

FIG. 14 is a flow chart showing a routine for 
calculation of a cylinder air filling amount Mci for a 
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second embodiment . 

FIG. 15 is a flow chart showing the routine for 
calculation of the value of a variable x. 

BEST MODE FOR WORKING THE INVENTION 
5 Below, embodiments of the present invention will be 

explained in detail with reference to the drawings. FIG. 

I show the case of application of the present invention 
to a four-stroke cylinder injection type spark plug 
ignition type internal combustion engine. Note that the 

10 present invention may also be applied to another spark 

plug ignition type internal combustion engine or a 
compression ignition type internal combustion engine as 
well . 

As shown in FIG. 1, in the present embodiment, an 
15 engine body 1 provided with, for example, eight cylinders 

is provided with a cylinder block 2, pistons 3 
reciprocating in the cylinder block 2, and a cylinder 
head 4 fixed on the cylinder block 2. Between each piston 
3 and the cylinder head 4, a combustion chamber 5 is 
20 formed. The cylinder head 4 is provided with, for each 

cylinder, an intake valve 6, intake port 7, exhaust valve 
8, and exhaust port 9. Further, as shown in FIG. 1, at 
the center of the corresponding part of the inside wall 
of the cylinder head 4, a spark plug 10 is provided. Near 
25 the inside wall of the cylinder head 4, a fuel injector 

II is provided. Further, the piston 3 is formed at its 
top surface with a cavity 12 extending from below the 
fuel injector 11 to below the spark plug 10. 

The intake port 7 of each cylinder is connected 
30 through an intake branch tube 13 to a surge tank 14. The 

surge tank 14 is connected through an intake pipe 15 to 
an air cleaner 16. Inside the intake pipe 15 is provided 
a throttle valve 18 driven by a step motor 17. Note that 
in this description, the part of the intake passage 
35 including the intake pipe 15 downstream of the throttle 

valve 18, surge tank 14, intake tube 13, and intake port 
7, that is, the part of the intake passage from the 



throttle valve 18 to the intake valve 6, will be called 
the "intake pipe part IM". On the other hand, the exhaust 
port 9 of each cylinder is connected through the exhaust 
branch tube and exhaust pipe 19 to an exhaust 
purification device 20 built in a catalytic converter 21. 
This catalytic converter 21 is communicated with the 
atmosphere through a muffler (not shown) . 

The intake valve 6 of each cylinder is driven to 
operate by an intake valve drive system 22. This intake 
valve drive system 22 is provided with a cam shaft and a 
switching mechanism for . selectively switching a 
rotational angle of the cam shaft to an advanced side and 
a retarded side with respect to the crank angle. When the 
rotational angle of the cam shaft is advanced, as shown 
by AD in FIG. 2, the opening timing VO and closing timing 
VC of the intake valve 6 are advanced, therefore the 
opening/closing timing is advanced. On the other hand, 
when the rotational angle of the cam shaft is retarded, 
as shown by RT in FIG. 2, the opening timing VO and 
closing timing VC of the intake valve 6 are retarded and, 
therefore the opening/closing timing is retarded. In this 
case, the amount of lift and working angle (opening time 
interval) of the intake valve 6 are maintained and the 
phase angle (opening timing) is changed. In the internal 
combustion engine shown in FIG. 1, the rotational angle 
of the cam shaft is switched to the advanced side or the 
retarded side in accordance with the engine operating 
state. Note that the present invention can be applied as 
well when the opening timing of the intake valve 6 is 
continuously changed or when the amount of lift or 
working angle are changed. 

Referring to FIG. 1, an electronic control unit 
(ECU) 31 is comprised of a digital computer provided with 
components connected to each other through a bi- 
directional bus 32 such as a RAM (random access memory) 
33, ROM (read only memory) 34, CPU (microprocessor) 35, 
input port 36, and output port 37. The intake pipe 15 



upstream of the throttle valve 18 is provided with an air 
flow meter 40 for detecting the flow rate of air (intake 
gas) passing through the intake pipe 15, Further, the 
surge tank 14 is provided with a pressure sensor 41 for 
detecting the pressure Pm of the air in the intake pipe 
part IM (below, referred to as the "intake pressure") . 
Further, the accelerator pedal 42 is connected to a load 
sensor 4 3 for generating an output voltage proportional 
to the amount of depression of the accelerator pedal 42, 
while the throttle valve 18 is provided with a throttle 
opening degree sensor (not shown) for detecting an 
opening degree of the throttle valve 18. The output 
signals of these sensors 40, 41, and 43 are input through 
corresponding AD converters 38 to an input port 36. 
Further, the input port 36 is connected to a crank angle 
sensor 4 4 generating an output pulse each time the crank 
shaft rotates by for example 30*". The CPU 35 uses the 
output pulse of the crank angle sensor 44 to calculate 
the engine speed. On the other hand, the output port 37 
is connected through corresponding drive circuits 39 to 
the spark plugs 10, fuel injectors 11, step motors 17, 
and intake valve drive system 22. These are controlled 
based on output from the electronic control unit 31. 

In the internal combustion engine of the present 
embodiment, the fuel injection amount (fuel injection 
time) TAUi of an i-th cylinder (i=l,2,..- ,8) is for 
example calculated based on the following equation (1). 

TAUi=TAUb-r|i- k (1) 
Here, TAUb indicates the basic fuel injection amount 
(basic fuel injection time) , r|i indicates the air amount 
variation correction coefficient of the i-th cylinder, 
and k indicates another correction coefficient. 

The basic fuel injection amount TAUb is the fuel 
injection amount required for making the air-fuel ratio 
match the target air-fuel ratio. This basic fuel 
injection amount TAUb is found in advance as a function 
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of the parameters relating to the engine operating 
conditions (for example, the engine load, the engine 
speed NE, etc., below referred to as the "operating 
parameters") and stored in the form of a map in the ROM 
5 34. Alternatively, it is calculated by an equation based 

on the operating parameters- Further, the correction 
coefficient k expresses the air-fuel ratio correction 
coefficient, acceleration increase correction 
coefficient, etc. combined together and is made 1.0 when 

10 no correction is needed. 

If the amount of air filled in an i-th cylinder when 
the intake stroke is completed is referred to as the 
"cylinder air filling amount Mci(g)", the air amount 
variation correction coefficient r|i is for compensating 

15 for variation in the cylinder air filling amount Mci 

among the cylinders. The air amount variation correction 
coefficient Tii of an i-th cylinder is for example 
calculated based on the following equation (2) . 
•ni=Mci/Mcave • • • (2) 

20 Here, Mcave indicates an average value of cylinder air 

filling amounts Mci (=EMci/8. Here, "8" expresses the 
number of cylinders) . 

For example, if deposits comprised mainly of carbon 
are formed on the inner circumference of the intake pipe 

25 part IM or on the outer circumference of the intake 

valves 6, since the amount of deposition of the deposits 
will differ for each cylinder, the cylinder air filling 
amount Mci is liable to vary among cylinders. Further, 
there is sometimes manufacturing error among cylinders in 

30 the volumes of the combustion chambers 5. In this case as 

well, the cylinder air filling amount Mci is liable to 
vary among the cylinders. In the case where the cylinder 
air filling amount Mci varies among cylinders, if fuel 
injection amounts are equal among all of the cylinders, 

35 the air-fuel ratios or the output torques will vary among 

cylinders. Therefore, in the present embodiment, the air 
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amount variation correction coefficient r\± is introduced 
to compensate for variations among cylinders in the 
cylinder air filling amount. 

Note that, considering the fact that the timing at 
5 which fuel is actually injected is ahead of the timing of 

calculation of a fuel injection amount TAUi by exactly a 
certain time, it is also possible to make the basic fuel 
injection amount TAUb in equation (1) a predictive value 
advanced by exactly a certain time from the timing of 
10 calculation of the fuel injection amount TAUi by equation 

(1) . 

Alternatively, the fuel injection amount TAUi of an 
i-th cylinder can be calculated based on the following 
equation (3) . 

15 TAUi=Mci- kAF- k (3) 

Here, kAF is a correction coefficient for making the air- 
fuel ratio match with a target air flow ratio. 

In this case as well, considering the fact that the 
timing when fuel is actually injected is ahead of the 

20 timing of calculation of the fuel injection amount TAU by 

exactly a certain time, it is also possible to make the 
cylinder air filling amount Mci in equation (3) a 
predictive value advanced by exactly a certain time from 
the timing of calculation of the fuel injection amount 

25 TAU. 

In this way, both when calculating a fuel injection 
amount TAUi based on equation (1) and when calculating it 
based on equation (3) , to make the air-fuel ratio match 
with the target air-fuel ratio for all of the cylinders 

30 so as to eliminate the variation of output torque among 

the cylinders, the cylinder air filling amount Mci has to 
be accurately found. 

In the present embodiment, a cylinder air filling 
amount Mci is calculated based on the drop of the intake 

35 pressure Pm caused by the intake stroke of the i-th 

cylinder, that is, the intake pressure drop APmdwni. 
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Next, the intake pressure drop APmdwni will be explained 
while referring to FIG. 3 to FIG. 5. 

FIG. 3 shows the intake pressure Pm detected by the 
pressure sensor 41 at for example constant time intervals 
5 over 720^ crank angle. The intake order in the internal 

combustion engine shown in FIG. 3 is #l-#8-#4-#3-#6-#5- 
#7-#2. In FIG. 3, OPi (i=l,2,... ,8) is the intake valve 
opening/closing timing of an i-th cylinder, while a 0° 
crank angle expresses the intake top dead center of the 

10 #1 cylinder #1. As will be understood from FIG. 3, when a 

certain cylinder starts an intake stroke, the rising 
intake pressure Pm will start to fall and therefore an 
upward peak will occur in the intake pressure Pm. The 
intake pressure Pm falls further, then again rises and 

15 therefore a downward peak will occur in the intake 

pressure Pm. In this way, the intake pressure Pm is 
alternately formed with upward peaks and downward peaks. 
FIG. 3 shows the upward peaks and the downward peaks 
caused in the intake pressure Pm due to the intake stroke 

20 of an i-th cylinder by UPi and DNi, respectively. 

As shown in FIG. 4, if the intake pressure Pm at an 
upward peak UPi is referred to as the maximum value 
Pmmaxi and the intake pressure Pm at a downward peak DNi 
is referred to as the minimum value Pmmini, the intake 

25 stroke of an i-th cylinder causes the intake pressure Pm 

to fall from the maximum value Pmmaxi to the minimum 
value Pmmini. Therefore, the intake pressure drop APmdwni 
in this case is expressed by the following equation (4) : 
APmdwni=Pmmaxi-Pmmini - • • (4) 

30 On the other hand, as shown in FIG. 4, if the intake 

valve 6 starts to open, the flow rate of air flowing out 
from the intake pipe part IM and sucked into a cylinder 
CYL, that is, the cylinder intake air flow rate mci 
(g/sec, see FIG. 5), starts to increase. Next, when the 

35 cylinder intake air flow rate mci becomes larger than the 

flow rate of air. passing through the throttle valve 18 
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and flowing into the intake pipe part IM, that is, the 
throttle valve air passage flow rate mt (g/sec, see FIG. 
5), the intake pressure Pm starts to drop. Next, when the 
cylinder intake air flow rate mci drops and becomes 
smaller than the throttle valve air passage flow rate mt, 
the intake pressure Pm starts to increase. 

That is, air flows into the intake pipe part IM 
through the throttle valve 18 by exactly the throttle 
valve air passage flow rate mt, and flows out from the 
intake pipe part IM through each intake valve 6 by 
exactly the cylinder intake air flow rate mci when air is 
taken into the i-th cylinder. In consideration of the 
above, the amount of outflow, that is, the cylinder 
intake air flow rate mci, temporarily exceeds the amount 
of inflow, that is, the throttle valve air passage flow 
rate mt, so the pressure inside the intake pipe part IM, 
that is, the intake pressure Pm, falls by exactly the 
intake pressure drop APmdwni. 

Now, the cylinder air filling amount Mci is the 
cylinder intake air flow rate mci integrated over time. 
Therefore, by ignoring the effect of overlap of the 
intake valve opening/closing timing OPi on the cylinder 
air filling amount Mci or air amount variation correction 
coefficient r|i (see FIG. 3), the cylinder air filling 
amount Mci can be expressed by the following equation 



Here, tmaxi indicates the timing when an upward peak 
occurs in the intake pressure Pm, that is, an upward peak 
generation timing, tmini indicates the timing when a 
downward peak occurs in the intake pressure Pm, that is, 
a downward peak generation timing, Atdwni indicates the 
time interval (sec) from the upward peak generation 



(5) . 




mt)dt + mt • 



Atdwni -h Atoc 



(5) 



2 
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timing tmaxi to the downward peak generation timing 
tmini, and Atoc indicates the intake valve opening time 
(sec) (see FIG. 4 ) . 

In equation (5) , the right side first term expresses 
5 the area of the part shown by Tl in FIG. 4 (below, 

referred to as the "region Tl"), that is, the area of the 
part surrounded by the cylinder intake air flow rate mci 
and the throttle valve air passage flow rate mt . The 
right side second term expresses the area of the part 
10 shown by T2 in FIG. 4 (below, referred to as the "region 

T2"), that is, the area of the part surrounded by the 
cylinder intake air flow rate mci, the throttle valve air 
passage flow rate mt, and the line mci=0, approximated by 
a trapezoid. 

15 As explained above, due to an intake stroke being 

performed, the cylinder intake air flow rate mci 
temporarily exceeds the throttle valve air passage flow 
rate mt. Therefore, the cylinder air filling amount Mci 
obtained by integrating the cylinder intake air flow rate 

20 mci over time exceeds the time integrated value of the 

throttle valve air passage flow rate mt . The region Tl, 
in this way, expresses the excess of the cylinder air 
filling amount Mci with respect to the time integrated 
value of the throttle valve air passage flow rate mt 

25 occurring due to an intake stroke being performed. 

Therefore, generally speaking, by dividing the 
cylinder air filling amount into a first amount of air 
expressed by the area of the region Tl (excess air 
amount) and a second amount of air expressed by the area 

30 of the region T2 (basic air amount) , making the first 

amount of air the excess of the cylinder air filling 
amount with respect to the throttle valve air passage 
amount occurring due to an intake stroke being performed, 
and totaling the first amount of air and the second 

35 amount of air for a cylinder, the cylinder air filling 

amount of the cylinder is calculated. 

On the other hand, the Law of the Conservation of 
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Mass for the intake pipe part IM can be expressed by the 
following equation (6) using the state equation for the 
air in the intake pipe part IM. 



dPm Ra - Tm 

5 = • (mt - mcx) • • • (6) 

dt Vm 

Here, Vm indicates the volume (m^) of the intake pipe 
portion IM, Ra indicates the gas constant divided by the 
average molecular weight of air (below, referred to as 
10 simply as the "gas constant"), and Tm indicates the 

temperature (K) of the air in the intake pipe part IM 
(see FIG- 5) . 

In the interval from the timing tmaxi to the timing 
tmini, the intake pressure Pm falls by exactly the intake 
15 pressure drop APmdwni. Therefore, by expressing 

Vm/(Ra-Tm) by the parameter Km and expressing the 
throttle valve air passage flow rate mt by the average 
value mtave, equation (5) can be rewritten using equation 
(6) to the following equation (7) . 



20 



A tdwni + A toe 
Mci = A Pmdwni • Km H- mtave • = • • • (7) 



This being the case, if detecting the intake 
pressure Pm by the pressure sensor 41 to calculate the 

25 intake pressure drop APmdwni, finding the above-mentioned 

parameter Km, detecting the throttle valve air passage 
flow rate mt by the air flow meter 40 to calculate the 
average value mtave, and detecting the timings 
tmaxi, tmini from the intake pressure Pm and throttle 

30 valve air passage flow rate average value mtave to 

calculate the time interval Atdwni (=tmini-tmaxi ) , 
equation (7) can be used to calculate a cylinder air 
filling amount Mci. Note that intake valve opening time 
Atoc is stored in advance in the ROM 34 . 
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However, as explained at the start, it is difficult 
to accurately find the intake pipe volume Vm and intake 
pipe temperature Tm. Therefore, in the present 
embodiment, the parameter Km is found without finding the 
5 intake pipe volume Vm and intake pipe temperature Tm. 

Next, a method for calculation of the parameter Km of the 
present embodiment will be explained while referring to 
FIG. 6A and FIG. 6B. 

In the present embodiment, note should be taken of 

10 the amount of air flowing into the intake pipe part IM 

and the amount of air flowing out of the intake pipe part 
IM in a set crank angle range set so that at least two 
cylinders for which the cylinder air filling amount Mci 
is to be calculated are included. 

15 FIG. 6A and FIG. 6B show the case where the 720*" 

crank angle range in which the intake strokes of all 
cylinders are included, for example, from intake top dead 
center of the #1 cylinder to the next intake top dead 
center of the #1 cylinder, is set as the set crank angle 

20 range. 

The total amount of air flowing into the intake pipe 
part IM within this 720° crank angle range is expressed as 
the area of the part shown by hatching in FIG. 6A, that 
is, the product of the throttle valve air passage flow 

25 rate average value mtave in this 720° crank angle range 

and the required time t72o required for the crank shaft to 
rotate by exactly 720° crank angle (mtave • t72o) - On the 
other hand, the total amount of air flowing out from the 
intake pipe part IM to the insides of the cylinders in 

30 this 720° crank angle range is expressed as the area of 

the part shown by hatching in FIG. 6B, that is, the total 
ZMci of the cylinder air filling amounts Mci. 

If the intake pressure Pm does not change much 
between the starting point and the end point in the 720° 

35 crank angle range, the total amount of air flowing into 

the intake pipe part IM in this 720° crank angle and the 
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total amount of air flowing out from the intake pipe part 
IM and filled in each cylinder in this 720® crank angle 
should be substantially equal to each other. Therefore, 
in this case, the following equation (8) stands. 

8 

mtave • t72o = ^ ^^i • • • (8) 

i = l 



Entering equation (7) at the right side of equation 
(8) and cleaning it up, the parameter Km can be expressed 
10 as in the following equation (9) . 

8 

2 t^2o - X tdwni + A toe) 

Km = mtave • ^ • • • (9) 

2 ^ A Pmdwni 

i = l 

That is, if calculating the throttle valve air 

15 passage flow rate average value mtave from the throttle 

valve air passage flow rate mt detected by the air flow 
meter 40, calculating the required time t72o from the 
output of the crank angle sensor 44, calculating the 
total value EAtdwni of the time interval Atdwni (see FIG. 

20 4) or the total value Z (Atdwni+Atoc) of the sum of the 

time interval Atdwni and the intake valve opening time 
Atoc (see FIG. 4), and calculating the total value 
ZAPmdwni of the intake pressure drop APmdwni, the 
parameter Km can be calculated. By doing this, the 

25 parameter Km can be simply found without finding the 

intake pipe volume Vm and intake pipe temperature Tm, 
therefore a cylinder air filling amount Mci can be simply 
and accurately found- 
As shown in FIG. 7A and FIG. 7B, for example, it is 

30 also possible to set a 360*^ crank angle range in which the 

intake strokes of four cylinders are included as the set 
crank angle range. In the example shown in FIG. 7A and 
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FIG. 7B, a first 360° crank angle range from intake top 
dead center of the #1 cylinder to intake top dead center 
of the #6 cylinder and a second 360° crank angle range 
from intake top dead center of the #6 cylinder to the 
5 next intake top dead center of the #1 cylinder are set. 

Regarding the first 360° crank angle range, the 
following equation (10) stands from the throttle valve 
air passage flow rate average value mtave at the first 
360° crank angle range, the required time tseo required for 

10 the crank shaft to rotate by exactly the first 360° crank 

angle range, and . the total SMcj (j=l,2,3,4) of the 
cylinder air filling amounts Mcj of the cylinders which 
perform an intake stroke in the first 360° crank angle 
range. Here, j indicates the order of the intake strokes. 

15 Similarly, regarding the second 360° crank angle range, 

the following equation (11) stands from the throttle 
valve air passage flow rate average value mtave' at the 
second 360° crank angle range, the required time t'360 
required for the crank shaft to rotate by exactly the 

20 second 360° crank angle range, and the total SMcj 

(j=5,6,7,8) of the cylinder air filling amounts Mcj of 
the cylinders performing an intake stroke in the second 
360° crank angle range. 



25 ^ Mcj = mtave • t^^^ • • • (10) 

3-1 



J] Mcj = mtave'-t'3go • • • (11) 

j=5 



Therefore, the parameter Km for the first crank 
30 angle range can be expressed by the following equation 

(12), while the parameter Km for the second crank angle 
range can be expressed by the following equation (13) . 
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Km = mtave • 



2 t. 



360 



- ^ (A tdwn j + A toe) 



(12) 



4 



2 ^ A Pmdwn j 
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Km = mtave' 



2 tV 



360 



(A tdwn j + A toe) 

j = 5 



(13) 



8 



2 ^ A Pmdwn j 



10 



15 



20 



25 



In this case, the cylinder air filling amounts Mcj 
(j=l,2,3,4) of the cylinders performing an intake stroke 
in the first crank angle range are calculated by equation 
(7) using the parameter Km calculated by equation (12), 
while the cylinder air filling amounts Mcj (j=5,6,7,8) of 
the cylinders performing an intake stroke in the second 
crank angle range are calculated by equation (7) using 
the parameter Km calculated by equation (13). 

Therefore, generally speaking, the set crank angle 
range is set so that the intake strokes of at least two 
cylinders for which the cylinder air filling amount is to 
be calculated are included, the total value ZAPmdwni of 
the intake pressure drops APmdwni of the cylinders 
performing intake strokes in this set crank angle range 
is calculated, and the above-mentioned first amount of 
air is calculated based on the intake pressure drops 
APmdwni and the intake pressure drop total value 
ZAPmdwni. Alternatively, it is also possible to view this 
as calculating the first amount of air based on the 
intake pressure drops APmdwni, intake pressure drop total 
value ZAPmdwni, throttle valve air passage flow rate mt 
or its average value mtave, required time required for 
the crank shaft to rotate by exactly the set crank angle 
range, time interval Atdwni from when the intake pressure 
Pm generates an upward peak UPi (see FIG. 4) to when it 
reaches a downward peak DNi or its total value ZAtdwni, 
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or intake valve opening time Atoc or its total value 
ZAtoc. 

However, for example, if the intake valve 
opening/closing timing is set to the retarded side RT 
5 (see FIG. 2), the intake valve closing timing VC will 

become the intake bottom dead center or later. In this 
case, even if the piston starts to rise, the intake valve 
6 is held in the open state, so the air sucked into the 
cylinder is liable to flow back into the intake pipe part 

10 IM. If this backflow occurs, as shown by X in FIG. 8, the 

cylinder intake air flow rate mci will temporarily become 
a negative value and the region T2 can no longer be 
approximated by a trapezoid- That is, when there is 
backflow of air from inside a cylinder to the intake pipe 

15 part IM at the end of the intake stroke, the cylinder air 

filling amount Mci cannot be accurately calculated from 
equation (7) . 

Therefore, in the present embodiment, when the 
intake valve opening/closing timing is set at the 

20 retarded side RT, the action of calculation of a cylinder 

air filling amount Mci by equation (7) is prohibited. In 
this case, the cylinder air filling amount Mci is not 
updated. The air amount variation correction coefficient 
T^i is calculated from the cylinder air filling amount Mci 

25 calculated at the previous calculation cycle. 

FIG. 9 shows the routine for calculation of the fuel 
injection amount TAUi of an i-th cylinder of the present 
embodiment. This routine is executed by interruption 
every preset crank angle. 

30 Referring to FIG. 9, at step 101, based on the 

engine load, the engine speed, etc. detected by the load 
sensor 43, the crank angle sensor 44, etc., the basic 
fuel injection amount TAUb is calculated. Next, at step 
102, the routine for calculation of the cylinder air 

35 filling amount Mci shown in FIG. 10 is executed. Due to 

this, the cylinder air filling amount Mci for a cylinder 
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is calculated. Next, at step 103, based on the cylinder 
air filling amount Mci calculated at step 102 and the 
average value Mcave of the cylinder air filling amount 
for all cylinders, equation (2) is used to calculate the 
5 air amount variation correction coefficient r|i of an i-th 

cylinder (i=l,2,... ,8). Next, at step 104, the 
correction coefficient k is calculated. Next, at step 
105, based on the basic fuel injection amount TAUb, air 
amount variation correction coefficient r\x, and 

10 correction coefficient k calculated at steps 101, 103, 

104, equation (1) is used to calculate the fuel injection 
amount TAUi. The fuel injector 11 of the i-th cylinder 
injects fuel by exactly the fuel injection amount TAUi. 
FIG. 10 shows the routine for calculation of the 

15 cylinder air filling amount Mci of an i-th cylinder in 

the present embodiment. 

Referring to FIG. 10^ at step 121, it is judged if 
the opening timing of the intake valve 6 is set at the 
advanced side AD (see FIG. 2). When the opening timing of 

20 the intake valve 6 is set at the advanced side AD, next 

the routine proceeds to step 122 where the throttle valve 
air passage flow rate average value mtave is calculated. 
Next, at step 123, the required time t72o is calculated. 
Next, at step 124^ the upward peak generation timings 

25 tmaxi and downward peak generation timings tmini for an 

i-th cylinder are detected (i=l,2,... ,8). Next, at step 
125, the time interval Atdwni of an i-th cylinder is 
calculated (Atdwni=tmini-tmaxi) . Next, at step 126, 
2 (Atdwni+Atoc) is calculated. Next, at step 127, the 

30 maximum values Pmmaxi and minimum values Pmmini of an i- 

th cylinder are detected. Next, at step 128, equation (4) 
is used to calculate the intake pressure drop APmdwni for 
the i-th cylinder. Next, at step 129, the intake pressure 
drop total value EAPmdwni is calculated. Next, at step 

35 130, equation (9) is used to calculate the parameter Km. 

Next, at step 131, equation (7) is used to calculate the 
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cylinder air filling amount Mci for an i-th cylinder. As 
opposed to this, when at step 121 the opening timing of 
the intake valve 6 is set at the retarded side RT, the 
processing cycle is ended. Therefore, calculation of the 
5 cylinder air filling amount Mci is prohibited. 

In the embodiments explained above, the region T2 
shown in FIG. 4 was approximated as a trapezoidal shape 
having a top side and a bottom side of Atdwni and Atoc. 
However, the region T2 may also be approximated as a 
10 rectangular shape with a side of for example Atdwni. In 

this case, the above-mentioned equations (7) and (9) 
become the following equations (14) and (15). 



15 



Mci = A Pmdwni • Km + mtave • A tdwni • * • (14) 



^720 " S ^ tdwni 

Km = mtave • q — ^ • • • (15) 

^ A Pmdwni 

iol 

Next, a second embodiment of the present invention 
will be explained. The equation (6) is modified as in the 
following equation (16) if expressing Vm/Ra as the 
20 parameter Km' . 

Km' dPm 

mt - mci = • • • • (16) 

Tm dt 

In the interval from the timing tmaxi to the timing 
25 tmini, the intake pressure Pm falls by exactly the intake 

pressure drop APmdwni, so the equation (5) can be 
rewritten using equation (16) to the following equation 
(17) . 

Km' Atdwni + Atoc 

30 Mci = APmdwni - + mt • • • • (17) 

Tm 2 
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This being the case, if detecting the intake 
pressure Pm by the pressure sensor 41 to calculate the 
intake pressure drop APmdwni, finding the above-mentioned 
parameter Km' , detecting the throttle valve air passage 
5 flow rate mt by the air flow meter 40 to calculate the 

average value mtave, and detecting the timings tmaxi, 
tmini from the intake pressure Pm to calculate the time 
interval Atdwni {=tmini-tmaxi ) , equation (17) can be used 
to calculate the cylinder air filling amount Mci. Note 

10 that in equation (17), the intake valve opening time Atoc 

is a value instructed from the ECU 31 to the intake valve 
drive system 22, therefore is the time when the intake 
valve 6 actually is open. 

However, as explained at the start, due to the 

15 backflow of the air sucked into a cylinder to the inside 

of the intake pipe or other factors, if calculating the 
cylinder air filling amount as explained above, the 
cylinder air filling amount ends up including an error. 
That is, the right side second term of equation (17) 

20 approximates the region T2 of FIG. 4 by a trapezoid. 

However, when backflow of air, etc. , occurs, the value 
calculated by approximation by the right side second term 
of equation (17) ends up greater than the region T2 by 
exactly the amount shown by the hatching in FIG. 11 and, 

25 as a result, the cylinder air filling amount is 

calculated greater and error ends up occurring. In other 
words, if making the intake valve opening time Atoc a 
value equal to the time when an intake valve 6 actually 
opens, the value calculated by approximation by the right 

30 side second term of equation (17) will end up including 

error. 

Therefore, in the present invention, by adjusting 
the intake valve opening time Atoc to a suitable value 
rather than making it the time during which an intake 
35 valve 6 actually opens, even if air flows back etc., the 

region T2 can be calculated with a high precision. Below, 
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while referring to FIG. 12A, FIG. 12B and FIG. 13, the 
method of calculation of the cylinder air filling amount 
in the present invention will be explained. 

FIG. 12A and FIG. 12B show the cylinder air filling 
5 flow rate mci and throttle valve air passage flow rate 

average value mtave for all cylinders in the crank angle 
720° from intake top dead center of the #1 cylinder to the 
next intake top dead center of the #1 cylinder. 

The total amount of air flowing into the intake pipe 

10 part IM during this crank angle 720"" is the area of the 

part shown by hatching in FIG. 12A, that is, is expressed 
by the product of the throttle valve air passage flow 
rate average value mtave during this crank angle 720® and 
the required time t72o required for the crank shaft to 

15 rotate by exactly the crank angle 720*^ (mtave • t72o) • On the 

other hand, the total amount of air flowing out from the 
intake pipe portion IM and being filled into the 
cylinders during this crank angle 720"^ is expressed as the 
area of the part shown by hatching in FIG. 12B, that is, 

20 the total ZMci of the cylinder air filling amount Mci. 

If the intake pressure Pm does not change much at 
all between the starting point and end point of the crank 
angle 720°, during this crank angle 720°, the total amount 
of air f lowing^ into the intake pipe part IM and the total 

25 amount of air flowing out from the intake pipe part IM 

and being filled into the cylinders should be 
substantially the same. Therefore, in this case, the 
following equation (18) stands. 

8 

30 mtave • t^^o = X! ' ' " ^-^^^ 



Further, if entering equation (17) at the right side 
of equation (18) and cleaning it up, it can be expressed 
as the following equation (19) . 
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g ^ ^ Atdwni + Atoc • 8 

mtave • t^^o = APmdwni • + mtave • — ' ' * (19) 

Tmave 2 

Here, Tmave indicates the air temperature average 

value in the intake pipe part IM during a crank angle 

720°. 

5 However, equation (19) may not stand in practice. 

This is because, as explained above, by making the intake 
valve opening time Atoc a value equal to the time for 
which an intake valve 6 actually opens, the right side 
second term of equation (17) causes an error to arise in 

10 the value calculated by approximation. 

Therefore, in the present invention, for the 
equation (19), the intake valve opening time Atoc is 
replaced by use of the variable x. In this case, the 
equation (19) is expressed as the following equation 

15 (20). 

8 

g ^ ^ Atdwni + X • 8 

mtave • t™ = Y APmdwni • + mtave • — * ' * (20) 

'^"^ f:i Tmave 2 

Further, by cleaning up equation (20) for the 
variable x, it can be expressed as the following equation 
20 (21). 

8 B 

Z Atdwni APmdwni 
Km' ^ 

^ _ ^720 _ 1=1 . 1^1 . . . (21) 

4 8 Tmave 4 • mtave 

The thus calculated variable x is a value 
25 corresponding to the intake valve opening time Atoc and a 

value determined when assuming that the total amount of 
air flowing into the intake pipe part IM during the crank 
angle 720° and the total amount of air flowing out from 
the intake pipe part IM and being filled in the cylinders 
30 during the crank angle 720° are equal (below, referred to 

as the "virtual intake valve opening time") . That is, the 
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virtual intake valve opening time x is a value set so 
that the area of the part surrounded by broken lines in 
FIG. 13 (that is, the trapezoidal part with a top side of 
Atdwni, a bottom side of the virtual intake valve opening 
5 time x, and a height of mtave) becomes equal to the area 

of the part surrounded by the cylinder intake air flow 
rate mci, the throttle valve air passage flow rate mtave, 
and the line mci=0 (region T2) . However, FIG, 13 shows 
the case for one cylinder, but in actuality, the virtual 

10 intake valve opening time x is set so that the total 

value for all cylinders of the areas of the parts 
surrounded by the broken lines becomes equal to the total 
value for all cylinders of the areas of the regions T2- 
On the other hand, equation (17) becomes the 

15 following equation (22) by using the variable x instead 

of Atoc: 

Km* ^ Atdwni + x 

Mci — APmdwni • + mt • • • • (22) 

Tm 2 

20 Further, by entering the value of the variable x 

calculated by equation (21) into equation (22), the 
cylinder air filling amount of each cylinder is 
accurately calculated. 

That is, according to the present invention, a 

25 virtual intake valve opening time whereby the average air 

flow rate to all cylinders becomes equal to the throttle 
air passage flow rate is calculated. By using the virtual 
intake valve opening time as the opening time of an 
intake valve to calculate the region T2, it is possible 

30 to accurately calculate the cylinder air filling amount 

to each cylinder. 

FIG. 14 shows a routine for calculation of the 
cylinder air filling amount Mci of an i-th cylinder 
according to a second embodiment. In the second 

35 embodiment, instead of the calculation routine shown in 

FIG. 10 of the first embodiment, the calculation routine 
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shown in FIG. 14 is used to calculate the cylinder air 
filling amount Mci of the i-th cylinder. 

Referring to FIG. 14, at step 141, a throttle valve 
air passage flow rate mt is detected from the output of 
5 the air flow meter 40 etc. Next, at step 142, from the 

output of the pressure sensor 41, the upward peak 
generation timing tmaxi and downward peak generation 
timing tmini of the intake pressure due to the intake 
valve 6 of an i-th cylinder opening are detected 

10 (i=l,2,... ,8). Next, at step 143, based on the peak 

generation timings tmaxi and tmini detected at step 142, 
the time interval Atdwni of an i-th cylinder is 
calculated (Atdwni=tmini-tmaxi) . Next, at step 144, the 
calculation routine of the virtual intake valve opening 

15 time x shown in FIG. 15 is used to obtain the calculated 

variable x. 

At step 145, from the output of the pressure sensor 
41, the maximum value Pmmaxi and minimum value Pmmini of 
intake pressure due to the intake valve 6 of an i-th 

20 cylinder opening are detected. Next, at step 146, based 

on the maximum value Pmmaxi and minimum value Pmmini 
detected at step 145, equation (4) is used to calculate 
the intake pressure drop APmdwni of the i-th cylinder. At 
step 147, based on the output of a temperature sensor 

25 (not shown) etc., the temperature Tm of the intake pipe 

part IM is detected- * Further, at step 148, based on the 
mt, Atdwni, x, APmdwni, and Tm calculated at steps 141, 
143, 144, 146 and 147, equation (22) is used to calculate 
the cylinder air filling amount Mci to a cylinder. The 

30 calculated cylinder air filling amount Mci to a cylinder 

is used for calculation of the fuel injection amount TAUi 
to the cylinder shown in FIG. 9. 

FIG. 15 shows the routine for calculation of the 
variable x according to an embodiment of the present 

35 invention. This calculation routine is performed each 

time the crank shaft rotates by 720"". 
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Referring to FIG. 15, at step 161, based on the 
output of the crank angle sensor 44 etc., the time t72o 
required for the crank shaft to rotate by 720° is 
calculated. Next, at step 162, from the output of the air 
5 flow meter 40 etc., the average value mtave of the 

throttle air passage flow rate while the crank shaft is 
rotating by 720° is calculated. Next, at step 163, the 
time interval Atdwni calculated at step 143 of FIG. 14 is 
totaled up for all of the cylinders to calculate ZAtdwni . 

10 At step 164, the intake pressure drop APmdwni calculated 

at step 146 of FIG. 14 is totaled up for all of the 
cylinders to calculate EAPmdwni. Next, at step 165, based 
on the output of the temperature sensor, the average 
value Tmave of the temperature in the intake pipe part IM 

15 is calculated. Next, at step 166, based on the t720/ 

mtave, SAtdwni, SAPmdwni, and Tmave calculated at steps 
161, 162, 163, 164, and 165, equation (21) is used to 
calculate the value of the variable x. 

Note that as explained above, equation (21) is 

20 conditional on the intake pressure Pm not changing much 

between the starting point and end point of the crank 
angle 720°, so it is preferable to calculate the cylinder 
air filling amount Mci only at the time of steady-state 
operation and to suspend the calculation of the cylinder 

25 air filling amount Mci at the time of a transitory 

operation where the intake pressure Pm easily fluctuates 
between the starting point and end point of the crank 
angle 720°. Here, the "steady--state operation" means for 
example the time of operation under an engine load or an 

30 engine speed which is substantially constant, while the 

"transitory operation" means for example the time of 
operation under an engine load or an engine speed which 
fluctuates . 

Further, in the above explanation, the present 
35 invention was applied to the case where the intake valve 

opening/closing timing was set to the retarded side, an 
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intake valve was open even after intake bottom dead 
center, and therefore the air sucked into a cylinder 
flowed back into the intake pipe. However, the present 
invention can be applied not only to this case, but also 
5 for example the case where the intake valve 

opening/closing timing is set to the advanced side, an 
intake valve opens from before intake top dead center, 
and due to this air does not flow into the intake pipe 
even while the intake valve is open. 
10 Note that while the present invention was explained 

in detail based on specific embodiments, a person skilled 
in the art could make various changes, modifications, 
etc. without departing from the claims and idea of the 
present invention . 



